Opto-Electrical Cooling of Polar Molecules 



ON 

o 
o 

(N 

> 

o 
O 



43 

i 



> 



^1" 
O 

o> 
o 



M. ZcppcnfcldQ M. Motsch, P.W.H. Pinkse, and G. Rempe 
Max- Planck- Institut fur Quantenoptik, Hans-Kopfermann-Str. 1, D-85748 Garching, Germany 

(Dated: October 7, 2009, PREPRINT v6) 

We present an opto-electrical cooling scheme for polar molecules based on a Sisyphus-type cooling 
cycle in suitably tailored electric trapping fields. Dissipation is provided by spontaneous vibrational 
decay in a closed level scheme found in symmetric-top rotors comprising six low-field-seeking rovibra- 
tional states. A generic trap design is presented. Suitable molecules are identified with vibrational 
decay rates on the order of 100 Hz. A simulation of the cooling process shows that the molecu- 
lar temperature can be reduced from 1 K to 1 mK in approximately 10 s. The molecules remain 
electrically trapped during this time, indicating that the ultracold regime can be reached in an 
experimentally feasible scheme. 

PACS numbers: 37.10.Mn, 33.80.-b 
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The ability to prepare samples of ultracold molecules 
opens up exciting new possibilities in physics and 
chemistry, including ultrahigh-precision molecular spec- 
troscopy and interferometry Q, investigations of 
anisotropic collisions in quantum-degenerate gases 
steering of chemical reactions tests of fundamental 
physics such as the search for the electron dipole mo- 
ment and novel approaches to quantum computing 
and quantum simulations @, 0, S]- Reaching ultracold 
temperatures through laser cooling has the great advan- 
tage that it does not lead to particle loss and that it 
is a single-particle process which does not require suit- 
able collision properties or high densities. However, laser 
cooling has so far only been demonstrated for atoms and 
ions with simple energy- level structures, whereas optical 
cooling of molecules has proven confoundingly difficult. 

Optical cooling of molecules requires a change in 
paradigm: In contrast to ultracold atoms, for which effi- 
cient cooling was realized early on || but trapping proved 
to be a challenge due to the shallow optical and magnetic 
potentials available, electric trapping of polar molecules 
is relatively easy [10( and has, in fact, been demonstrated 
for molecules without any cooling Optical cooling 

of molecules could therefore start with trapped molecules 
and exploit the tremendous (~ 1 K) energy-level shifts 
producible by laboratory electric fields, circumventing 
the usual requirements of standard laser cooling such as 
highly closed transitions, fast decay rates, and significant 
photon momentum transfer. 

Making use of the aforementioned paradigm shift we 
here present a cooling scheme for molecules which is con- 
ceivable with present technology. Specifically, we replace 
photon recoil by an electric-field interaction energy as 
the means to remove energy from a molecular ensemble 
in a configuration reminiscent of Sisyphus cooling and 
single-photon cooling H, H, 14|. Spontaneous emission 
of photons serves only to remove entropy. As a result, 
the number of scattered photons required to achieve sub- 
stantial cooling is dramatically reduced. Slowly decaying 



vibrationally excited states, generally offering stricter se- 
lection rules than electronic transitions, can therefore be 
used for the spontaneous decay. 
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FIG. 1: (color online). Energy-level and state-transition di- 
agram for the cooling scheme. A molecule in the strong lfs 
state |s) diffuses from the low- field region to the high-field 
region (1) where it is driven to the weak lfs state \w) (2). 
After moving back to the low-field region (3), the molecule 
is driven to the excited state |e) (4) from which it decays 
spontaneously to \s) (5). The irreversible spontaneous decay 
makes this cycling process unidirectional. 

Our cooling scheme is shown in Fig. [TJ Two neighbor- 
ing regions in space, each with a constant but different 
electric field, are realized by a suitable arrangement of 
electrodes. These electrodes also provide a high-electric- 
field enclosure around these regions to ensure trapping of 
molecules in low- field-seeking (lfs) states. Fig. [2] shows a 
possible design for the electrodes. Trapped molecules ex- 
perience a potential step when moving from one region to 
the other. The magnitude of this potential step depends 
on the average orientation of the electric dipole moment 
of a molecule with respect to the electric field and may 
vary significantly for different molecular states. For one 
strong and one weak lfs molecular state we obtain a po- 
tential as a function of position as depicted by the curves 
|s) and \w) in Fig. [TJ 

Suppose a molecule possesses an excited state |e) which 
decays into the states \w) and \s). We induce transitions 
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FIG. 2: Design of an electric trap for the cooling scheme. Re- 
gions of tunable homogeneous fields are achieved using paral- 
lel capacitor plates (a). Collisions with the plate surface are 
eliminated by alternatingly-charged microstructured surface 
electrodes [15[ (b). Transverse confinement is achieved by a 
high voltage electrode between the plates around the perime- 
ter of the trap (a). By interrupting the perimeter electrode, 
an electric quadrupole guide can be connected to the trap for 
the injection and extraction of molecules [llj (c). 

between \w) and |e) in the low-field region of the trap and 
transitions between |s) and \w) in the high- field region. 
Doing so creates a unidirectional cycling process. During 
the cycle, the molecule loses a kinetic energy correspond- 
ing to the difference between the potential steps of the 
strong and the weak lfs state, leading to overall cooling. 

The main advantage of this cooling scheme is the large 
amount of kinetic energy which can be removed from a 
molecule for each spontaneously emitted photon. For a 
representative dipole moment d e ; of 1 Debye [D] , oriented 
in an electric field E of lOOkV/cm, one obtains an inter- 
action energy of d e ; • E = |/cb x 1.61 K. Starting with an 
ensemble of molecules with a translational temperature 
below 1 K, this in principle allows the removal of all of 
a molecule's kinetic energy in a single step. In practice, 
however, more than one spontaneous decay is necessary 
to cool a molecule: When the fields are kept constant, a 
molecule will generally end up in state \s) in the low- field 
region with insufficient energy to move back to the high- 
field region but with at least the amount of energy ob- 
tained when moving from the high- to the low-field region 
in the state \w). Further cooling to lower temperatures 
therefore requires the height of the electric-field step to 
be slowly ramped down, allowing the cooling cycle to re- 
peat. Nonetheless, a few dozen spontaneous decays are 
more than enough to cool a molecule to below a mK. 

Due to the small number of spontaneous photon emis- 
sions, the requirements imposed on the emission process 
are much less stringent than for standard laser cooling. 
Not only is the branching ratio for decay from the excited 
state to desired and undesired states much less critical, 
but the rate at which such transitions occur may also be 
much lower. As a result, the use of vibrational transitions 
for the spontaneous decay process is possible. 

The advantage of vibrational compared to electronic 
excitations is that except in the case of strong reso- 
nances with other vibrational states, a molecule with one 



molecule 


del 


/vib 


7 


molecule 


del 


/vib 


7 




[D] 


[cm" 1 ] 


[Hz] 




[D] 


[cm" 1 ] 


[Hz] 


CFH 3 


1.85 


2964 


37 


CF3CI 


0.50 


1105 


73 


CF 3 H 


1.65 


3036 


65 


CF3B1- 


0.65 


1089 


74 


CH 3 CCH 


0.78 


3334 


87 


CF3I 


0.92 


1080 


61 


CF 3 CCH 


2.36 


3327 


79 


BH3CO 


1.80 


2217 


274 


N(CH 3 ) 3 
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TABLE I: An overview of symmetric-top molecules with 
strong parallel vibrational transitions with permanent dipole 
moment d e i [lq |. transition frequency / v ib [13], and sponta- 
neous decay rate 7. The italicized values were obtained using 
the quantum-chemistry package Gaussian fl8|. We have suc- 
cessfully produced a cold sample of each of the molecules on 
the left using our quadrupole guide [19[|. Note that the large 
hyperfine splitting in CF3CI, CFsBr, and CF3I complicates 
the straightforward application of the present scheme to these 
molecules. 



quantum of excitation in a single vibrational mode will 
decay primarily back to the vibrational ground state. 
Additionally, compared to the deep ultraviolet wave- 
lengths required to excite electronic states of most simple 
chemically-stable molecules, many molecules have strong 
vibrational transitions in the wavelength range 3—10 /jm. 
The coverage of this wavelength range by tunable narrow- 
band light sources has been significantly improved in re- 
cent years by the commercial availability of quantum- 
cascade lasers and optical parametric oscillators, in ad- 
dition to, e.g., lead-salt lasers. 

Beyond the closed vibrational transition, the rota- 
tional transitions must be considered. The excited vi- 
brational state must not only decay to a manageable 
number of rotational states, but each of these states 
must be lfs so that the molecule remains trapped. Disre- 
garding linear molecules due to their generally weaker 
quadratic Stark interactions, symmetric-top molecules 
have the most stringent selection rules for dipole tran- 
sitions. Describing the rotational states of a symmetric- 
top molecule by the quantum numbers for the total an- 
gular momentum J, the angular momentum about the 
molecule's symmetry axis K and the angular momen- 
tum about a lab-fixed axis M , the selection rules for 
a parallel transition are AJ = 0, ±1, AK — and 
AM = 0,±1 [2(J. Furthermore, the lowest-order Stark 
interaction is E staTk = -E • d el = -\E\\d e i\ : j^^. Ob- 
serving that J > 0, |A'| < J and \M\ < J [2£J], we see 
that an excited state with \K\ = J > 2 and M = —K 
best satisfies the conditions stated above. Such a state 
may decay into a total of only five rotational states, all of 
which are lfs. These can be repumped using additional 
lasers or microwave fields. Note that the condition of few 
decay channels to purely lfs states can also be satisfied 
for linear molecules, e.g., in a £ electronic state using a 
vibrationally excited state with M — 0, J > 3 [20| . 

Use of vibrational excitations for opto-electrical cool- 
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FIG. 3: (color online). Transition scheme used to simulate 
cooling of CF3H. Rotational states are labeled using the no- 
tation \J, K, M). MW and ir denotes the induced microwave 
and infrared transitions. V denotes spontaneous decay to the 
five states with vibrational excitation v= 0. The energy lev- 
els are obtained by diagonalizing the rigid-rotor Hamiltonian 
for CF3H at field strengths of 5kV/cm and 20kV/cm for the 
low-field and high-field region, respectively. 



ing requires molecules with a sufficiently fast vibra- 
tional spontaneous decay rate. Table U lists promis- 
ing symmetric-top molecules. Although a decay rate of 
~ 100 Hz is glacial relative to decay rates used for laser 
cooling of atoms, it is adequate considering the small 
number of decays needed for the scheme presented here. 
Nonetheless, the spontaneous decay rate raises the ques- 
tion how fast opto-electrical cooling progresses. This is 
studied by numerically solving rate equations for cool- 
ing of CF3H. The rate equations and their derivation 
are included in the appendix. The first excitation of 
the C-H stretch mode at 3036 cm" 1 in the rotational 
state J = K = — M = 2 is used as the excited state. 
The fact that this state spontaneously decays to five 
v= states (v being the vibrational quantum num- 
ber) necessitates a somewhat more complicated transi- 
tion scheme than the one shown in Fig. [1] Specifically, 
we simulate cooling using the transition scheme shown 
in Fig. [3] The IR transition as well as each of the mi- 
crowave transitions are driven with a rate of 10 kHz. As- 
suming a Stark-broadening of 10 MHz, this would require 
narrow-linewidth sources with an intensity on the order 
of lmW/cm 2 for all the transitions involved. Sponta- 
neous decay from the excited state is modeled using a rate 
of 65.2 Hz, partitioned among the states with v= based 
on rigid-rotor dipole-transition matrix elements [20| . The 
volumes of both trap regions are set to 100mm 3 , con- 
nected by an area of 10 mm 2 . 

At time t = 0, molecules are distributed among the 
states v= in both trap regions with a v 2 dv velocity dis- 
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FIG. 4: (color online). Velocity distribution in the weak-field 
region of the trap after cooling for a time t. Velocities are 



converted to temperatures according to 
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itoring the population in the excited state during the cooling 
process shows that on average a molecule spontaneously de- 
cays 4.7x during the first second, 17. Ox during the first 5 
seconds, and 9.0 x during the next 5 seconds. 



tribution up to a cut-off velocity of 11.7m/s. This is the 
maximal trappable velocity of the involved states due to 
higher-order Stark shifts. The electric-field-strcngth dif- 
ference between the two trap regions is a free parameter 
which is varied as a function of time to be proportional to 
the 80th percentile of the kinetic energy of the molecules. 
The potential-energy step for each of the molecular states 
is modeled using the first-order Stark shift -Estark- 

The rate of the cooling process is influenced by three 
effects. Most significantly, the rate coefficients indicate 
the time in which 1 — 1/e of molecules perform some pro- 
cess, whereas the time in which 99% of molecules per- 
form this process takes significantly longer. Ramping 
down the electric-field step too rapidly therefore causes 
the final energy of most molecules to substantially exceed 
the field-step energy so that efficient cooling is no longer 
possible. Secondly, the fraction of energy removed during 
each cooling cycle is below unity. Reducing the temper- 
ature by, e.g., a factor of ten requires several cooling cy- 
cles. Finally, spontaneous decay to the states |2, 2,— 1), 
|3, 2, —1), and |3, 2, —2) in the low-field region of the trap 
has no net effect, reducing the effective decay rate. 

The velocity distribution of the molecules in the low- 
field region of the trap for various times after cooling 
commences is shown in Fig. SJ As can be seen, significant 
cooling occurs in under a second. Note that the cooling 
rate decreases significantly as time progresses. For high 
temperatures, the cooling rate is limited by the decay 
rate of the vibrationally excited state, allowing the tem- 
perature to decrease exponentially with time. For low 
temperatures, the cooling rate is limited by the time it 
takes for the molecules to move between the two regions 
of the trap, with the cooling rate proportional to the 
velocity of the molecules. Therefore, at very low temper- 
ature the cooling process is no longer efficient, and the 
molecules must either be moved to a smaller trap or a 
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different cooling scheme must be applied. 

The elementary description of opto-electrical cooling so 
far glosses over several issues which must be addressed to 
ensure the experimental viability of the method. In par- 
ticular, achieving required trapping times, sufficient mix- 
ing of the individual velocity components, and validity of 
approximate selection rules are now discussed. 

In addition to collisions with the background gas, Ma- 
jorana flips and rovibrational heating by thermal black- 
body radiation are the identified loss channels for po- 
lar molecules stored in electric traps [2l], [52]. Al- 
though rotational heating is a problem for extremely 
light molecules [22j ] and vibrational heating for heavy 
molecules, neither is the case for the molecules considered 
in table [H For example, the heating rates to the lowest 
vibrational modes never exceed a few mHz at 300 K for 
CF3CCH, the heaviest molecule in table HI 

Majorana flips are expected to have been a problem 
in past trap designs with a near-zero electric field in the 
central trap region [2l[. However, the trap in Fig. [2] is 
specifically designed to allow a homogeneous offset field 
throughout the vast majority of the trap volume. Fur- 
thermore, field zeros near the edges of the trap can be 
reduced to singular points through clever electrode de- 
sign, which essentially eliminates Majorana flips. 

Opto-electrical cooling only removes energy from a sin- 
gle component of the velocity vector, making sufficient 
mixing of the velocity components a necessity. Electric- 
field inhomogeneitics near the microstructured plate sur- 
face allow such mixing on a sufficiently short timescale. 
This is demonstrated by trajectory simulations discussed 
in the appendix. 

The zero-field rigid-rotor harmonic-oscillator selection 
rules used so far imply a closed six-level system for 
opto-electrical cooling. These selection rules are modi- 
fied in several ways for real molecules. Transitions with 
AK 7^ and decay to other excited vibrational states 
are generally possible for symmetric top molecules via 
resonances between near-degenerate vibrationally excited 
states. Due to the few spontaneous emissions needed, 
such couplings will at most cause problems for individ- 
ual molecule species. 

For non-zero values of the electric field, J ceases to be 
a good rotational quantum number and spontaneous de- 
cay with |A J\ > 2 becomes possible. The resulting con- 
sequences were checked by diagonalizing the rigid-rotor 
Hamiltonian for non-zero electric fields using molecular 
constants of CF3H and calculating dipole transition ma- 
trix elements between the new eigenstates. Although the 
partitioning of spontaneous decay from the state v= I , 
|2, 2, — 2) to the five states v=0, J = 2 and 3 is signifi- 
cantly changed already at electric fields of ~ 50kV/cm, 
the spontaneous decay to states with J > 4 remains be- 
low 1% for fields up to lOOkV/cm. This effect on opto- 
electrical cooling is therefore negligible. 

Achieving a temperature below 1 mK through opto- 



electrical cooling would allow other cooling schemes, re- 
quiring longer interaction times or higher phase-space 
density, to be implemented. Specifically, opto-electrical 
cooling can easily be extended to an accumulation 
scheme, for example to load molecules into a tightly con- 
fining optical dipole trap. The low temperatures and high 
densities thus achieved create extremely favorable start- 
ing conditions for a number of further cooling schemes 
such as evaporative cooling, cavity cooling, or sympa- 
thetic cooling with ultracold atoms. 
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DERIVATION OF THE RATE EQUATIONS 

The rate equations used to simulate cooling of CF3H 
can be derived as follows. The ensemble of molecules 
in the trap is represented by the number of molecules 
Pa (v,t)dv with velocity between v and v + dv in the 
molecular state a in part i = 1, 2 of the trap. Here, i = 1 
denotes the low-field region and i = 2 denotes the high- 
field region of the trap. 

In this description of the molecular ensemble, we ignore 
the position of the molecules within each trap region as 
well as the direction of the velocity vector v with v = |v|, 
effectively assuming instantaneous spatial redistribution 
of the molecules within each trap region as well as instan- 
taneous redistribution of the direction of v. Whereas the 
assumed instantaneous spatial redistribution should at 
most slightly affect the validity of the results of the sim- 
ulation since the molecule's thermal motion will rapidly 
redistribute the molecules through the trap, the assumed 
mixing of the components of v is less obvious. In fact, in 
the case of a potential energy in the trap which is com- 
pletely separable in Cartesian coordinates, no mixing of 
the velocity components would occur at all. 

To address this question, molecule trajectory simula- 
tions in a trap based on the design in Fig. [2] of the paper 
were performed and the temporal correlation of the mag- 
nitude of the individual components of v was calculated. 
The calculated correlation for a particle with a velocity 
of lOm/s in the homogeneous-field region of the trap is 
shown in Fig. [5l Particularly for large r, the correlation 
function for v x and v y can be accurately reproduced by 
assuming a ~ 20 % probability for the velocities to com- 
pletely mix for each collision with the microstructured 
plate surface. This demonstrates that significant mixing 
of the velocity components occurs due to the field inho- 
mogencities near the plate surfaces. Due to the trans- 
lational symmetry of the microstructures along the z- 
direction, this mixing does not include the z-component 
of the velocity, leading to significantly slower mixing for 
this velocity component. However, by arranging the mi- 
crostructure plates such that the structure on the top 
and bottom plate is rotated by 90 with respect to the 
other, all three velocity components mix on the shorter 
timescale. The resulting mixing is sufficiently strong that 
a major impact on the cooling rate does not occur. 

To derive the rate equations for the molecule distribu- 
tions pi 1 •* (v, t), two processes must be taken into account. 
First, transitions between the various internal molecular 
states may occur in both regions of the trap. This is 
modeled by a fixed fraction of molecules in a given state 
switching to a different state per unit time interval, 



d_ 

dt 



P «M) = £(c^W) 



(i) 




FIG. 5: Temporal correlation ((\vi(t)\ - \vi(t + t)\) 2 } of the 
individual components of the velocity for a particle in a trap 
based on the design in Fig.[2]of the paper. The inset shows the 
orientation of the velocity components relative to the trap. 



with appropriate rate coefficients c^\,. 

Second, molecules may diffuse between the two trap 
regions. For those molecules in region i with a veloc- 
ity component v x perpendicular to the interface between 
the trap regions, a fraction of v x A/Vi of the molecules at- 
tempt to enter the other trap region per unit time. Here, 
A is the surface area of the interface between the trap 
regions and Vi is the volume of trap region i. For those 
molecules in region 2 of the trap, all molecules attempt- 
ing to enter region 1 succeed, whereas for molecules in 
region 1, only those molecules with v x > 2AE a /m suc- 
ceed, where AE a is the potential energy difference be- 
tween the two trap regions for molecules in state a and 
m is the molecular mass. Finally we need to average over 
the possible values of v x . For an isotropic velocity distri- 
bution in three dimensions, a single velocity component 
is evenly distributed so that the fraction of molecules 
leaving region 2 per unit time is 



A v A 

— / dV-r.V-r— = - — , 



2v 



V 2 4 V 2 



(2) 



and the fraction of molecules leaving region 1 per unit 
time is 



1 

2^ 



.4 



dv x v T — 



v 2 „, A 



Vi Av Vi 



(3) 



Here, v a ,2 = \A' 2 — 2AE a /m is equal to the velocity a 
molecule in region 1 will have once it has reached region 
2. Note that the integrals are normalized by dividing by 
2v since v x can be both positive and negative but only 
molecules with (in this case) positive v x can enter the 
opposite trap region. 
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Every molecule which leaves one trap region must enter the opposite region. Nonetheless, a complication arises 
since the velocity of the molecules changes when they move between the regions. Molecules which enter region 1 in 
state a with velocity in the range v to v + dv must have had a velocity in the range v a ^ to v a .2 + vdv/ v a ^ in region 2. 

As a result, the number of molecules in region 2 which could potentially increase the number of molecules p^\v, t) dv 

(2) 

in region 1 with velocity in the range v to v + dv by switching regions is equal to p a (v a ,2,t)vdv/v aj 2- Combining 
this with Eqs. ([2|) and ((3]) one obtains a diffusion rate 



d_ 

dt 



¥| x ^Pa\v a , 2 ,t)- ^p^(v,t), v > ^2AE a /m 
£Hv,t)={ . (4) 



diii 



0, V < y/2AE a /m 



Note that the v in Eq. ([2]) is the velocity in region 2 and is therefore replaced by v a ^ in Eq. 0. For v 2 < 2AE a /m, 
the diffusion rate is zero due to the potential energy step. 

For molecules entering region 2 we introduce v a .\ = ^Jv 2 + 2AE a /m and the derivation of the diffusion rate is 
entirely analogous. One obtains as the final result, 



TZ> = x(&] a p£Hv,t)-c™,p£\v,t)), v < ^AEJ^ (5) 



These rate equations are used for the cooling simulations in the main text. 



